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Abstract 

We consider systems of higher spin gauge fields that are described by a free 
field Lagrangian and one interaction of arbitrary order N that is local and sat- 
isfies abelian gauge invariance. Such "solitary" interactions are derived from 
Noether potentials solving the Noether equations. They are constructed using 
a free conformal field theory carried by the same flat space as the higher spin 
fields. In this field theory we consider N-loop functions of conserved, confor- 
mally covariant currents, they are UV divergent. The residue of the first order 
pole in the dimensional regularisation approach to the iV-loop function is a lo- 
cal differential operator and is free of anomalies, so that current conservation 
and conformal covariance is maintained. Applying this operator to the higher 
spin fields, the Noether potential results. We study the cases N = 2, N = 3 
and N = 4. We argue that our N = 3 vertex for any number of derivatives A 
is identical with the known cubic interaction. 
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1 Introduction 



In the investigations of higher spin gauge field theories that started with the sem- 
inal work of Fronsdal [1], the trilinear or cubic interactions of such fields played a 
central role. For small spins, i.e. two and three, Berends, Burgers and Van Dam [2] 
constructed such interactions. In their work also the idea of finite or infinite series of 
iV-linear interaction Lagrangians was put forward and this was connected with the 
concept of higher order gauge transformations, meaning that a k-th order contribu- 
tion to a gauge transformation depends on k higher spin gauge fields locally, and 
on the gauge parameter e linearly. Till today no such series has been constructed 
beyond N = 3. We consider here only the gauge transformations for the case k = 0, 
called "abelian". 

After this classical work a long period of time was filled with investigations on the 
existence of general cubic interactions [3] and group theoretical analyses of criteria 
on the three spins which allow the interactions to exist [I]. After these results 
became known, the first explicit examples of such interactions were produced, the 
examples, however, belong to classes with special properties, e.g. conserved vector 
currents were constructed from two fields (such as the energy-momentum tensor) 
and these were coupled to the third field [5j [6]. Finally explicit expressions for 
general cubic interaction Lagrangians were found [H El E], and it was shown that 
the same expressions appear in the analysis of singular string models |10j . Recent 
results are summarized in [TT] . [T2] and [T5] . 

Both the cubic interactions (N = 3) and the free Lagrangian (N = 2) play an 
important role in the present investigation. Their explicit expressions are used to 
check (and understand) the results of the quantum field theory approach to any 
degree N interaction Lagrangian with abelian gauge invariance obtained here. We 
prefer the notation hS s \x; a) for the higher spin gauge field on the flat space with 
coordinates x where, with an arbitrary tangent vector a at x, the totally symmetric 
rank s tensor field is homogeneous of degree s in a. This notation allows to 
express all local operations on these tensor fields by differential operators in x and 
a. 

The results on cubic interactions as presented in El E] can be summarized as 
follows. The interaction Lagrangian consists of a sum of expressions jC^'f^f/i^], < 
i,j < 3, i + j < 3 each of which is trilinear in the fields 



h^\x; a) 



Dh (si) (x;a 




and has the form 



3-j-i 



(1.2) 
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where V fc denotes k spacetime derivatives distributed in any way among the three 
types of fields, and the total number of derivatives k + i is (Metsaev's inequality jl]) 

3 

Sj — 2 min{gj} < k + i (1.3) 

i=l 

This total number is fixed for a Noether potential and denoted A. Those expressions 
not containing deDonder or trace fields are denoted "leading" terms. From the 
derivation of the cubic interactions we learned several facts: 

1. For i = only the two cases C^'^ 3 ,j G {0, 3} are nonzero; 

2. Consequently in deDonder gauge, where Dh^ = 0, h and h decouple as 
dynamical variables; 

3. "Proposition 1" is valid: All £^'£3 can be reconstructed from the leading terms, 
they are unique. 

Proposition 1 follows from the construction technique used in [8]. Of course gener- 
alizations for arbitrary iV would be of great interest. 

The strategy used in J7J El E] has the following crucial aspects: 

1. In the coordinate representation the cubic interaction density is independent 
of the space dimension D. D is a "hidden" parameter; 

2. The existence and value of Metsaev's minimal differential operator degree A m j n 
jl] is taken as granted. For each A > A min a unique interaction Lagrangian 
(up to an overall normalization factor that may be chosen to depend on D of 
course) is obtained. 

The approach discussed in this work is based on a free scalar conformal quantum 
field theory carried by the same space as the higher spin fields. We construct con- 
served conformal currents each of an even spin s from this scalar field, and evaluate 
the iV-loop function of such currents. This approach has analogous aspects as the 
ones presented above from El E]. We have a nontrivial D dependence of the 
loop graph function in momentum space, and its UV divergence is of course also 
D dependent, so that dimensional regularization can be applied. In the resulting Y 
functions D-dependence of the pole positions is obvious, and the sum of the residues 
which gives the singular part of the loop function we are interested in, is D depen- 
dent. But the final integration over the loop momentum Q cancels the essential part 
of this dependence, namely up to its nonessential normalization factors (see (2.41), 
(2.42) and the comment following these equations). This allows us then to turn to 
part 2 of the strategy. We choose now A > A min freely as is done in [7J El [9] . In 
each case this has a consequence on the value of D. But part 1 of the strategy is 
interpreted as the permission to continue constant functions of D "analytically" in 
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D, allowing any integral value for D, even negative ones. We end up with the con- 
clusion that for any A we find a unique cubic vertex function and this is supposed 
to be identical with that in [3 El IS]- However, only the case A = A min is studied 
in detail in Section 5. For all terms calculated in Section 5 this identity between 
the results of both approaches has been verified. The construction gives a further 
result, namely that the vertex functions for any N > 3 are conformally invariant, 
when conformal covariance of the higher spin functions is appropriately defined. For 
the case N = 2, A = 2 the free Lagrangian is reproduced modulo trace terms (in 
the deDonder terms (1.1) the traces O a h^ must also be neglected). For N = 4 we 
have also constructed the leading terms. 

We define Noether potentials that are functionals of the higher spin gauge fields 
of the integral form 

QVp[hM] = f dxJ*\h^}(x) (1.4) 



and are defined to be scalar polynomials of degree N in derivatives of K higher spin 
fields and to obey 

1. locality, 

2. invariance with respect to abelian gauge transformations of all their functional 
variables off shell. 

Here we consider a system of K such fields 

{h^(z;a)}f =1 (1.5) 
The abelian gauge transformations are defined by 

5 h {s) (x; a) = (aV)e (s " 1} (x; a), □ a e (s - 1) (x; a) = (1.6) 
and abelian gauge invariance of the Noether potential is expressed by 

(Vd a ) 5 -^^(x-a) = (1.7) 

for all s and off-shell. These are the Noether equations. There are obviously N of 
these, one for each field contained in Q^. An interaction Lagrangian of degree N 
denoted Cjq[h] has the same locality properties but abelian gauge invariance is now 
fulfilled if the traceless part of 

(V(9 a ) ^ - | n-shell (1-8) 

vanishes. 

The Noetherian solving (1.7) has a polynomial form and is easily decomposed 
into the sum of two terms. The first vanishes on shell. We denote it the "kernel". 
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The remainder is the candidate for the iV-th order Lagrangian £jv[/i ] if iV > 3. 
However, it may be useful to apply field redefinitions to this to obtain a standard 
form. Such field redefinitions are not investigated in this article. 

In Section 2 we introduce the free scalar conformal quantum field theory which 
is used as a tool for our construction. The conserved currents of this one-field 
model are of particular interest. They have all even spin. Among these there are 
conformal currents and nonconformal derivatives of these. We argue that only the 
conformal currents play a role in our construction. Moreover, in order to treat a 
higher spin field theory of Fronsdal's type, these fields are doubly traceless. The 
corresponding conformal currents should then be doubly traceless themselves, too, 
to eliminate unwanted extra interactions of the trace fields. Such currents are given 
in Section 2, eqs. (2.12), (2.14)-(2.16) and contain in the general case one free 
parameter multiplying a traceless current inside the doubly traceless current. They 
can and must be employed in our algorithm as long as we consider the local covariant 
distribution defining the vertex. However, if we contract this vertex distribution with 
external fields that are doubly traceless, the double and higher trace subtractions 
of the currents drop out, but the free parameter multiplying a simple trace current 
remains. With a trick we can show that the free parameter can also be transferred 
to the higher spin field where it multiplies a simple trace term of the field. So the 
QFT calculations can always be performed with the naive currents (2.2). The free 
parameter ("B-parameter") enters the vertex function again by means of the higher 
spin fields to which a simple trace term with this parameter is added at the end (see 
(2.18)). The result of the 2- vertex has identical leading terms with the free Fronsdal 
Lagrangian. On the other hand, what concerns the check of our algorithm for the 
N = 3 vertex with the known result pH EJ |9] , we study here only the case A min and 
compare the leading terms and some divergence terms in both cases and find that 
they agree. Thus for the check on this level one can set the trace fields in the vertices 
for N = 2 and N = 3 equal zero by hand, and use the naive currents (2.2). We shall 
follow this strategy for N = 4, too. Since Proposition 1 fixes unique expressions 
for trace and deDonder terms from the leading terms, (and suggesting that it holds 
also for N = A) there cannot remain a free parameter. We will argue later that 
the following Proposition 2 is valid: The iV-point vertex functions and consequently 
the interaction Lagrangians for N > 3 can be calculated by our algorithm with the 
naive currents (2.2). 

Loop graphs for N conserved conformal currents of arbitrary spin need a renor- 
malization which is performed by dimensional regularization. The sum over the 
residues of all first order poles of the T-functions at positions n G No (there is only 
a finite number of them) gives a local differential operator which produces the iV'th 
order Noether potential. These residues being free of anomalies, current conserva- 
tion and conformal covariance remain valid, and consequently the Noether equation 
is satisfied. The fact that the parameter D, the space-time dimension, drops out and 
A is fixed instead, implies that the counting label for the UV poles n is no longer 
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limited by an expression (2.35) containing D, but instead simply (see (2.37)) by 

< n < — (1.9) 

In Section 3 the parameter integrals that define the residues are further evaluated, 
leading to formulae to which we refer in the sequel. 

Section 4 is devoted to the N = 2 free Lagrangian, leading to a first check of our 
method. Section 5 is devoted to the vertex for N = 3. As we emphasized, this is in 
principle a check of our quantum field theoretic ansatz. In Section 6 we study the 
case N = 4. If N > 3 typical complications arise due to the fact that iV momentum 
vectors submitted to momentum conservation possess bilinear forms that cannot 
be expressed linearly by momentum squares (Laplacians). This problem is dealt 
with in this article in a specific way, some elementary and technical relations are 
presented in an Appendix. The freedom of choosing this approach is equivalent to 
the possibility to act with partial integrations on the interaction Lagrangian density. 

2 The tools of the construction 

2.1 The conserved currents 

Given a free quantum field theory of a scalar massless field 

□$(x) = (2.1) 

Spin s currents are obtained for even s from this scalar field by an operator product 
expansion 

j (s) (x; a) = (a, Vi - V 2 ) s $(a;i)$(a; 2 ) | x x = x 2 = x (2.2) 
and their divergence is 

(Vi + V 2 ; d a )(a, V x - V 2 ) s $0n)$(:r 2 ) = 
s((a, Vi - V 2 ) s - 1 (n 1 - HkWxjQfa) = (2.3) 

since the free field $(x) satisfies the massless Klein-Gordon equation. For a related 
use made of massless free scalar fields and their conserved currents see [13] . 

The conformal currents of spin s are therefore all conserved. However, it is less 
well known that there are more conserved currents jn^ of spin s which are derivatives 
of conformal operators but not all conformal themselves 

ji s \x; a) = [(a 2 V 2 - (a, V) 2 ]™ x j^ 2n \x; a); n G {0, 1, 2, ..s/2} (2.4) 

Conservation of these currents follows from 

[(V,d a ),(a 2 V 2 -(a,V) 2 )] = (2.5) 
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It follows then that the 2-point function of any pair of such different currents vanishes 
since 

< | / s - 2ni) (:zi; a 1 )j ( - s - 2n2 \x 2 ; a 2 ) | >= if m ^ n 2 (2.6) 

because of the inequality of the two conformal representations (spins and conformal 
dimensions are different). 

We define the contraction of two symmetric tensors of equal spins by the star 
operation 

j^(x; a) * a h^(x; a) = (s\)- l j^{x; d a )h^(x; a) (2.7) 
and a generating function by the vacuum expectation value < | \P | > of 

#[yM ; fcM] = exp J d D xj {s) (x; a) * a h^ s \x; a) (2.8) 

where the vacuum belongs to the free conformal field theory of the scalar field 
operators $(a;). This generating function can be expanded in powers N of h starting 
with the trivial 1 for iV = and then a second order (N = 2) term. Each power 
iV yields after regularization and extraction of the singular part the iV'th order 
Noetherian. Since the exponent of the generating function is second order in the 
free field $(x), the vacuum expectation value of \1/ is feasible. But it is unclear of 
how to regularize and extract the singular terms without using a series expansion of 
the exponential function in advance. The vacuum expectation value of the current 
j' 8 ' vanishes by definition (normal ordering). 

If instead of the conformal current j( s > one would use the general conserved 
current 

s/2 

jW(x;a) = £4jW(x;a) (2.9) 

n=0 

one would destroy conformal covariance on the one hand, and on the other hand 
introduce factors of powers of 

aV-(a,p) 2 (2.10) 

where p is an external momentum, Laplacians and double divergences acting on 
the field h^ s \ In the interaction Lagrangians there is no place for such terms (as 
known from the case N = 3, [7]), they should be eliminated by field redefinitions if 
possible. In the case N = 2 on the other hand the number of derivatives is restricted 
to A = 2. The orthogonality relation (2.6) requires that at least two such factors 
(2.10) appear if A n ^ for n > 0, implying four and more derivatives. These are 
unwanted and therefore excluded. So it seems that the conformal covariant current 
is the only reasonable candidate for the generating function. 

Since the current j^(x; a) in (2.8) is conformally covariant (spin s, dimension D+ 
s — 2) we can define conformal covariance properties of h^ s \x; a) by the requirement 
of conformal invariance of the generating function, i.e. the representation of 
is dual to that of jW. Then any unrenormalized iV-loop function produced by the 
generating function is also conformally invariant. 
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The conserved currents (2.2) possess nonvanishing traces denoted by 

4 s \x;a)=/ s \x;a) 

j[ s ~ 2) (x; a) = — — ^—-j (s) (x; a) 
s(s - 1) 

^ 4> ^ a )= s ( S -l)(°-2)( S -3) ^ a) ete ' (2 - U ) 

Since the operators D a and (V, d a ) commute, conservation and vanishing of multiple 
traces are compatible requirements. Therefore starting from the conserved currents 
j( s \x;a) we can construct multiple traceless conserved currents by the expansion 
("subtraction of traces", B = 1) 

J«(x;o) = J2 B n ] (a 2 r.Jt 2n) (x;a) (2.12) 

n>0 

Requiring the single trace of J^(x; a) to vanish leads to0 

B^(sin 9 le) = (-D" 22 „„,|jt;-21„ - < 2 ' 13 > 

If the double trace of J^ s \x; a) is required to vanish, we use B[ s) = B^ (double) as 
a free parameter and expand 

(double) = + B[ s) p(*\ a { s) = p[ s) = 1, a[ s) = ^ = (2.14) 

Solving the recursions in this case yields 

„W = (_l)«-i !ft n>2 (2.15) 

V ' 2 2n n(n - 2)\[s + fj, - 3]„' ~ V ; 

^ 1 ij 2 2 «- 2 (n-l)![ S + / i-4]„_ 1 ' - 1 l ^ ibj 

If in (2.7) we replace p s \x\ a) by either of these choices j( s >(x; a), then only that part 
of the field h^ s '(x;a) which is traceless (respectively is doubly traceless) is coupled 
to the current. By comparison of (2.16) with (2.13) we find that 

^ = B^? (single) (2.17) 

belonging to a simply traceless current. Though (1.8) expresses that the divergence 
of the traceless part of the current must vanish, this requirement is trivial if the 
total current is conserved. Thus b[ s ^ remains a free parameter indeed. Of course 
it is also trivial that any doubly traceless current remains doubly traceless if an 
arbitrary simply traceless current is added. 



1 We use the notation [x] n — 11^=0 ( x ~ ^) an< ^ A* = D/2, D the space dimension 
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Contracting the current j( s \x;a) with the doubly traceless higher spin field 
h^(x;a) by the * a -operation (2.7) gives 



{j (s) {x . a)+ l a 2 {Daj) (s-2) {x . a)} ^ h {s) {x . a) = 
[S\2 

n(<) 

j (s \x; a) * a {h (s \x; a) + -^- a 2 (D a h) (s ~ 2 \x; a)} (2.18) 

[ S \2 

Consequently in all our formulas for the vertex functions we can evaluate N-point 
loop functions as vacuum expectation values of the currents j^ s \x;a) (2.2) and 
multiply these loop functions not with but with + a 2 Trh^ s ~ 2 \ In this 
article we calculate all terms including trace and deDonder terms only for N = 2, 
the leading plus some divergence terms for N = 3 and only the leading terms for 
N = 4. In this context we refer again to the Proposition 1 formulated in Section 
I. Its validity for N > 4 has not been proved yet, however, but it could be true in 
general. 

2.2 The N-point loop functions 

In the Hilbert space of the free field $(x) over the D dimensional flat space we 
consider the connected iV-point function 

N 

<0\]Jj isi) (x i ;a i )\0> conn (2.19) 
i=i 

which consists of 

±(N-1)\ (2.20) 
contractions (by Wick's theorem) e.g. of the type 

N D 
HiaiWiY^ia^W^Y^dxi-x^) 2 )-^, fi = - (2.21) 

i=i 

We restrict our discussion to this example. All other examples would be mathemat- 
ically equivalent. 

Going to the momentum representation 

/ Y[ dxj exp{iXjPj}, / Y[ dxj exp{i( Xj - x j+1 )q jJ+1 } 

J 3=1 J 3=1 

Vi -> Pi = q i: i + i - q t - lt i (2.22) 
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and neglecting constant and only D-dependent factors, and postponing a factor 
we have the integral 

N N N 

n^M+i(?M+i) _1 Yl( a mj+iY j ~ rj H(a k q k -i,k) rk 
i=i i=i fc=i 

AT 

j J^'i/'/ - + (2.23) 

Here labels 1 < % < N are used cyclically, so that i + 1 for £ = N means 1 and k — 1 
for k = 1 means N. Moreover from the Fourier integrals imaginary units % result 
which multiply to (— 1) R i s , R = J2i r «> $ = J2i s i that cancel the r-dependent signs 
from all expansion coefficients = (-l) r '{ s ;) in (2.2). 

These delta-functions can be taken into account by solving for the (fr^+i in terms 
of the Pi and a new momentum Q 

N 

= Q + (pi +P2 + ••• + Pi) = Q + w u qN,i = Q, ^2pi = w N = (2.24) 



i=i 



The resulting integral is 

N 

/ dQ[Q 2 (Q + Wl ) 2 (Q + w 2 ) 2 ...(Q + wn^) 2 ]- 1 Y[(oi(Q + w t )) s ^(a t (Q + w^))^ 
J i=i 

(2.25) 

The square bracket in the integral can be represented by 

/ ]T exp{-[aiQ 2 + a 2 (Q + wi) 2 + a 3 (Q + w 2 ) 2 + ... + <tjv(Q + w N ^) 2 ]} x ... 
Jo 4 

(2.26) 

We introduce 

N 

1=1 

so that the exponent can be rewritten as 

{-S[(Q + r 2Wl + t 3 w 2 + ... + r^jv.i) 2 + Z(t, w)}} (2.28) 

Now denote 

N-l 

Q' = Q + J2 Ti + lWi = Q + *( r > w ) ( 2 - 29 ) 
i=i 

and introduce the integration variable Q' instead of Q. 
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This must be done also in the gradient terms after multiplication with | 
(now the sign factor from (2.23) is taken into account) and summing 



N , s. N 

J2H(r) ( a '(G + ^) St ~ n «Q + ^-i)' n = i[«2Q + Wi + (2.30) 

{n} i=i ^ ri ' i=i 

and then inserting Q' 

2Q + Wi + Wi = 2Q' + 2Ri, R i ( y T,w) = ^( y w i + w i - 1 )-x(r,w) (2.31) 

(a,, Q' + ifc)* = ( St ) M) mi ^iRi) Si - mi (2.32) 

In the sequel we skip the factor 2 s , S = J2 i Sj. 

Now we rescale Q' by Q' = Yr^Q" and obtain after Q" integration 

poo 

/ dSI! JV - 1 -3(^+ D )ex P {-I!Z} 
= T(N --(M + D))Z~ N+ ^ M+D) (2.33) 

N 

M = ^m,, (2.34) 

i=l 

In order to obtain a local vertex, 

n = I(M + D)-iV (2.35) 

must be natural nonnegative, because this and only this guarantees that this result 
is polynomial in {wi}. Therefore we must have an ultraviolet divergence from the 
gamma function. In its argument we replace D — > D — 2e leading to a pole behaviour 

1 (— l) n 

F(N — -(M + D)) = - — p-e _1 + 0(1) (2.36) 

and we have to collect all terms with first order pole in e and sum over their residues. 
Since Z is of second order in the momenta (derivatives), 2n is the number of deriva- 
tives contributed by Z n . The power of momenta contributed by the second factors on 
the r.h.s. of (2.32) is S — M. Thus in order to obtain a fixed number A = S — M + 2n 
of derivatives of the whole expression we must allow an n-dependence of M: replace 
M by M n = M + 2n, and correspondingly we use A n = S — M n = A — In. From 
the last equation follows that n is bounded: 

n < ^A (2.37) 
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so that only a finite number of terms survive. Using only the singular e _1 terms 
guarantees that conservation of the free field currents from which we started is not 
violated by anomalies. Moreover the conformal covariance of these terms is also 
maintained. This is the source from which the conformal invariance of the vertex 
functions follow. There is another consequence: since n is an integer M + D must 
be even, i.e. M is restricted respectively to even (odd) values if D is even (odd). 
The integration over q = Q" 

N 

/ dqe^{-q 2 }\{M m \ (2-38) 
J i=i 

is somewhat involved. The trick consists in introducing 

N / M \ N 

5>M] M = £( rUrMr (2.39) 

tT tr \mi,m 2 ,m 3 ...m N J f = \ 

and integrating 

dqexp{-q 2 }(uq) M = {u 2 )^ M / dQ D ^ cos M 9 / dpp M + D ~ x exp{-p 2 } 

N 

Y^&i, \q\=p (2-40) 



N 

u 

i=l 



If M is odd, the l.h.s. vanishes since the integrand changes sign under q — > —q. 
Thus from now on M is assumed to be even, implying also that D must always be 
even (see (2.35)). The remaining spherical integral and the integral over p are 



/ 



Jo 



dQn^ cos M 9 = 2n^ F ^ M+1) l (2.41) 

r(i(M + £>)) v ; 

1 1 

dpp M+D ^ exp{-p 2 } = ^(i(M + D)) (2.42) 



From (2.41), (2.42) we recognize that an explicit dependence on D (apart from 
trivial factors) drops out of our formulae after multiplication of (2.41) with (2.42). 
We can freely continue our results analytically in D, e.g. for iV = 3, Si = s 2 = s 3 = 
4, A = 4, M = 8 we may continue to the value D = —2. 

Now we expand, neglecting in this first study systematically all the trace terms 
proportional to a 2 (as was discussed in Section 1) and defining n^j = for i > j, we 
get from (2.40) 

( M 2)M/2 = 2 m/ 2 ( M/2 '\ J] $f^ +n ^ + trace terms (2.43) 
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so that by comparison with (2.39) 



mi = y^X n hk + n k,i) (2-44) 
k 

and the coefficient of Yli j( a i a j) ni ' j * s 

The remaining integrations are over r iy i G {2,3, ...N}. They are all elementary. 
We call n the "order" of the terms. For a first illustration of our algorithm we 
study mainly the "leading terms" for fixed n, this means that we neglect trace 
and divergence (or deDonder) terms. For arbitray order n the factor Z n yields 
polynomials of degree n in the differential operators PiPj which for N = 3 (and only 
in this case) can be expressed by Laplacians pf. We emphasize that the leading 
terms are explicitly independent of D, this follows from the multiplication of the 
two integrals (2.41), (2.42). The traceless and the doubly traceless currents ((2.12)- 
(2.16)) depend on D explicitly by the parameter jj, = y. 



3 Evaluating the residues 

For the leading terms we perform all integrations over the parameters r. We start 
from (2.31), (2.32) where the dependence on the r is in the Ri and from Z, the 
quadratic complement in the exponent (2.28) 

N N N 

Z = ^Tiw1_i - ^TiWi-if = ^ Ti^Wi-uWj-!) (3.1) 
i=2 i=2 i,j=2 

where 

Ti,j = Ti(5ij - Tj) (3.2) 
We use that the integral over the Tj can be extended to include a variable t\ by 



/ 



N N 

n^(i-$>,) (3.3) 

i=l j=l 



which allows us to express Z for iV = 3 in the beautiful symmetric form 

Z = 7iT 2 (pi) 2 + t 2 t 3 { P2 ) 2 + r 3 r 1 (p 3 ) 2 (3.4) 

This implies that the terms of order n contribute an n'th order polynomial of the 
Laplacians. For iV > 3 contractions of gradients (pi,Pj),i ^ j, cannot be eliminated. 
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We emphasize that Euler integrals of the type 

/ (TT ^ r , ^ - y Ti) = v h (3.5) 

that are easily derived by induction, are the optimal tools for the integration over 
Tj. Integrals related with these are 

/N N yrN , 

iUdTjT? 6{t 3 )}6{1 - J» = ( ^n+ nV (3-6) 

Now we turn to the Ri terms (2.30)-(2.32). We expand the expression (2.32) in 
trinomial form 



x(a,g') mi (« 4 , + Wi-JWoix) 1 * (3.7) 

where 

rji = Si - rrii - h (3.8) 
Besides Z n only the last factor in (3.7) has to be integrated over the Tj. 

3.1 The case n = for general N 

We consider now the simplest case, namely order zero 

n = (3.9) 
Then from the last factor in (3.7) and using (2.29) we get the integral 

N N , N N , v JV-1 

'H*rj)H.( a *) Vi = /(II^IBE J\ nfw^r-} (s.io) 

3=2 i=l J i=2 i=l {tt^} Vl^^i/ a=1 

where the exponents 7i ita form partitions of rji, without order and with zeros included. 
We introduce the shorthand ( a with the properties 

N N-l 



^ ^ ~^i,a Con ^ ^ ~^i,a Vii 
i=l a=l 

JV-1 AT 

C« = - mi - ki) = S - M - K, K = Y,ki (3.11) 



q=1 i=l 
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From (3.6) and (3.10) we get by performing the integration 

Multiplying with the trinomial coefficient from (3.7) we can write the coefficients as 

Si \f Vi \ ( ^ , (313) 



ITLi 



Since M is fixed for fixed n (here zero), we obtain contributions where M/2 is 
decomposed into a partition of numbers {n i; j} forming a triangular matrix: n it j = 
for i > j (see (2.43)). For any given set {n^} the summation parameters {rrii} (or 
{me}) are fixed by (2.44) 



partitions {mj} VL MJ/ VL J/ i,J 

nJV-1 > , N , 

_ g=i w TTf-iv»M s * 

rr ^ (S-M-K + N-l)\ l\ { ) U*<>{^} 

partitions {iri ya },k£ t=l 



N-l 



(hat^we + we-!)))"' \[{a,w a y^ (3.14) 



a=l 



The sum over a of the partitions 7Tj Q gives by (3.11) the variable r^. The £ Q are by 
(3.11) considered as functions of these partitions in the same way. 

The last two factors of (3.14) represent a polynomial of w, respectively a differ- 
ential operator of degree 

A = S-M, S = J2 S * ( 3 - 15 ) 



3.2 The case n = 1 for general N 

In the case n = 1 we assume that the number of derivatives is the same, but two of 
them are reserved now for Z. Therefore we have to replace 

A = S-M by A 1 = A-2 = S-M 1 = S-M-2 (3.16) 

and in (3.14) M by Mi = M + 2 everywhere and then insert the factor Z 
t-tAt-i >. | n^-Vi 



(S-M-K + N-iy. J (S-M-K + N-l)\ 

N-l N-l 

x{{S - M - K + N - 2) + 1)^ - E(C/J + 1H] 2 } (3-17) 

0=1 0=1 
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Instead of (3.14) we obtain then 




M/2 + 1\ (M + 2 



) 



i 





0=1 

JV-1 N , x 

-£«, + iwn(-.)«^; { J 




(3.18) 



7=1 

In general n extends from zero to y (see (1.9), (2.37)). 



4 The free Lagrangian and the 2-vertex 

The general treatment of the Noether potential for arbitrary N in Section 3 can be 
applied to the case N = 2, too. In this case the two higher spin fields must be identi- 
cal by QFT. It sounds reasonable to use the double traceless currents (2. 12), (2.14)- 
(2.16) including also the B-parameter, to check the trace terms of the 2-vertex 
against the trace terms of the free Lagrangian. We obtain immediately 



So we have obtained one differential operator in the tensorial spaces of a±, a 2 applied 
to one vacuum expectation value. If this 2-point function is contracted with two 
equal higher spin fields of spin s, which are double traceless, then in (4.1) all terms 
in the differential operator containing powers (a\a%) n , n > 2 drop out. On the 
other hand a[ s ^ vanishes by (2.14). Therefore the whole differential operator in (4.1) 
reduces to the unit operator. If we allow next for the B-parameter, we use (2.18) to 
resume it in the higher spin field (see below). 

Since A = 2 admits only the two poles at n = and n — 1, we can simply carry 
over the results of Section 3. The main generalization we must admit concerns the 
matrix n it j of exponents of (a^, a,j) which is upper triangular still but now with a 



< | J {s \x l -a l )J {s \x 2 -a 2 ) | > 



= J2 a n )2 (*lai) n < | Jt 2n) (x 1 ;a 1 )Jt 2n \x 2 ;a 2 ) | > 



n=0 




(4.1) 
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nonvanishing diagonal for the traces, for N = 2 

n id = if % > j 

7^1 = 2711,1+711,2, 7722 = 2712,2 + 711,2 (4.2) 

But we have to modify the eqs (3.14), (3.18) in the powers of 2: Replace in (3.14) 
2 M / 2 and in (3.18) 2 M / 2+1 by the same expression 2 E »<; ni + 

Then in (3.14) we consider the last two factors. There is only one momentum 
variable contained in them 

w 1 =p 1 = -p 2 = v (4-3) 

and these factors take the form 

2-^+ k ^(a 1 ,p) kl+7Vl - 1 {a 2 ,p) k2+7T2 ' 1 = 2-( kl+k2) (a 1 ,p) s ~ mi (a 2 ,p) s - m2 (4.4) 



4.1 The case n = 

From A = 2 follows M = 2s — A = 2s — 2. For the exponents {fiij} the following 
list of possible values exists 

7ll,l 77 2 ,2 77i,2 777i m 2 h + 7Ti,i k 2 + 7T 2 ,1 

(1) s-1 s-1 s-1 1 1 

(2) 1 s-2 s s- 2 2 (4.5) 

(3) 1 s-2 s-2 s 2 

(4) 1 1 s - 3 s - 1 s - 1 1 1 

For each of these four cases (l)-(4) we have the following possible values of k± : 7Ti,i, k 2 , 7r 2 ,i 

ki k 2 7Ti,i 7T 2 ,l 

(1),(4) 110 

(1),(4) 10 1 

(1),(4) 110 

(1),(4) 11 

(2) 2 (4.6) 

(2) 10 1 

(2) 2 

(3) 2 
(3) 10 1 
(3) 2 

These we have to insert into the modified formula (3.14) 

r(i(2s-i)) y 2 n ^- k 4 5-1 V 2a_2 V 

2 uu \ni,i,n 1>2 ,n 2:2 J \m 1 ,m 2 J 



Til, 1,^1,2,712, 2, ^1,^2,^1,1,^2,1 

xf-1 



7Tl 1+ ^ 2 X (^1,1 + ^2,1) 



(3 - fcx - k 2 )\ \m 1 , k 1: n 1A J \m 2 , k 2 , n 2 ,i, 
x (ajy^ (al) n ^ (ai, a 2 ) ni - 2 K p)*"™ 1 (a 2 , p)^ 2 (4.7) 
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The result consists of three terms 



^^{4s(ai,p)(a 2 ,p)(ai,a 2 ) s 1 

+s(s - l)[(a 2 ,p) 2 a\ + (a 1} p) 2 al] (a l} a 2 ) s ~ 2 
+s{s - l)(s - 2)(ai,p) (a 2 ,p) (a 1? a 2 ) s " 3 a? a 2 } (4.8) 

4.2 The case n = 1 

The derivatives are now all in 

K) 2 = P 2 (4.9) 

so that 

fcl + 71-1,1 = h + 7T 2 ,i = (4.10) 

It follows that only two possibilities appear 



lin n 22 fi 12 f s ) 

' > > Vni.l, 712,2, "1,2/ 

(a) s 1 

(/3) 1 1 s-2 s(s-l) 



(4.11) 



which yield the expression 

" s! V ^ p 2 [4(ai, a 2 ) s + s(s - l)( 0l , a 2 ) s - 2 a\ a 2 ] (4.12) 



4! 2 

Thus we obtain for the sum of the terms (4.8) and (4.12) 

^^-{[4s(a 1 ,a 2 ) s ~ 1 + s(s - l)(s - 2)(a 1 ,a 2 ) s ~ 3 a 2 1 a 2 l }(a 1 ,p)(a 2 ,p) 

+s(s - ^[(ai.pfa 2 , + (a 2 ,p) 2 a 2 ] (a± : a 2 ) s ~ 2 
-p 2 [4(ai, a 2 ) s + s{s - l)( Ql , a 2 ) s ~ 2 a\a 2 2 }} (4.13) 

The free Lagrangian is on the other hand 

£o[h (s) ] =~h^(x;a) * a ^ s \x;a) 
+ 8s{s 1 _ 1) D a h^(x- a) * a □ a ^W(x; a) (4.14) 

where the Fronsdal term is defined as 

J*"\x; a) = nh^ s \x; a) - {aV){Vd a )h^ s \x; a) + UaV) 2 U a h {s \x; a) (4.15) 
U a ^ s \x- a) = 2[UU a - (Vd a ) 2 ]h^(x; a) (4.16) 
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In (4.13) we must replace a, by d ai and apply the resulting differential operator to 
h^ s \x; ai)h^(x', a 2 ). 

First we notice that the free parameter b[ s ^ enhances the problem with the trace 
terms if it does not vanish. Namely with the trace field h( s ~ 2 ^ = [s(s — l)]~ 1 D a /i <s ' ) 
for a doubly traceless field we get 

+ -^a*n a h^] = [1 + 2 ^*ip 4 W ' } (4-17) 

and we obtain a coefficient depending on the space dimension D. Such coefficients 
are unknown for N = 2 or N = 3. 

However, if we neglect the trace terms hS B ~ 2 \ the Fronsdal field (4.15) reduces 

to 

J* s) (x; a) \ U s-2) =0 = Dh (s \x; a) - (aV)(V<9 a )/i (s) (x; a) (4.18) 



and (4.13) simplifies to 

(s!) 3 v^F 



6-2 

which is, up to the normalization, the same as (4.14). 



h {s) (x;a) * a J r (x;a) \- h (s )=0 (4.19) 



5 The cubic Noetherian 
5.1 The leading terms 

In [7j the leading terms of the cubic Lagrangian (those which remain when the trace 
and deDonder fields vanish) are presented in the cyclic form 

= Y, C ni%Z j d Zl dz 2 dz 3 5(z 3 - Zl )5(z 2 - z x ) 

x (d a o b ) Q ^ (d b d c ) Q ^ (d c d a ) Q ^ (d a v 2 r (d b v 3 r 

x(«9 c Vi) n3 /i (si) (,2i; a)h {s2 \z 2 ; b)h^\z 3 ; c) (5.1) 
The total order of space derivatives is 

A = ni +n 2 + n 3 (5.2) 
and the equations for the tensorial orders are 

n\ + Qi, 2 + Qz,\ = si, n 2 + Q 2)3 + Q 1>2 = s 2 , n 3 + Q 3A + Q 2t3 = s 3 (5.3) 
These equations are solved by 

Q2,3 = nt - v u Q 3A =n 2 - u 2 , Q\,i = n 3 - u 3 , ra» > v { > (5.4) 
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implying 

1 3 
Vi = Si - -(S - A), S = ^2 s i ( 5 - 5 ) 

i=t 

Since from [4] we know that the minimal order A is 

A m j„ — S — 2s 3 in the case s\ > s 2 > s 3 (5.6) 

we obtain that in this case 

v\ = st - s 3 , v 2 = s 2 - s 3 , v 3 = (5.7) 

From now on we consider only the case of minimal A min (5.6). Finally we give the 
coefficient 



5.1.1 The case n=0 

We concentrate for simplicity on the example A min . In the cubic case the exponents 
m, and n^j form triplets and these can be mapped one-to-one onto each other and 
on the variables Qij and rij. These relations are 



"l2 


= Ql2 


= "3 


- v z 


= s 3 


- m 3 




= g 3 i 


= n 2 


- v 2 


= S 3 


- m 2 


"23 


= Q23 


= "i 


- V\ 


= s 3 


— rrii 



i>i = -(A + Si — Sj — Sk), i,j, k different (5.9) 

so that the trinomial coefficient in (5.8) agrees with the trinomial factor in (2.43). 
These relations represent the first step towards a reconstruction of the structure of 
the cubic Lagrangian. 

Now we consider the w-polynomials contained in (3.14), (3.18) leading to the 
differential operators. We neglect all divergence expressions of the type 

(a iPi ) = (5.10) 

together with all trace terms proportional to any a 2 

~(ae,w £ + wi-t) = {ae,wi-i) (5.11) 
and in detail we get by inserting Wi from (2.9) 

(ai,toi) = 0, {a 2 ,wt) = -(a 2 ,p 3 ), (a 3 , wt) = (a 3 ,pi) 
(at,w 2 ) = {at,p 2 ), {a 2 ,w 2 ) = -(a 2 ,p 3 ), (a 3 ,w 2 ) = 

(5.12) 
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Consequently in (3.14) we obtain 

ki = h = 0, 7Ti,i = 7r 3)2 = (5.13) 

implying 

7Ti,2 = 7?i = Si - mi, 7T 2 ,i + 7T 2 ,2 = ??2 = S 2 ~ m 2 - fc 2 , ^3,1 = ^3 = S 3 - m 3 (5.14) 

Thus the w-polynomial comes out as 

(ai,p 2 ) Sl - mi («2,p 3 ) S2 - m2 (a3,Pi) S3 - m3 (5.15) 

which leads to the differential operator expected from the balance equations and 
the cyclic ordering. Thus there remains only to check the numerical factor. The 
coefficients in front are 

33 )( V (5.16) 

"1,2, "2,3, "1,3/ \ m l, m 2, m 3 J 

The desired final trinomial coefficient is the first factor of (5.16), which was stated 
already after (5.9), meaning that the second trinomial coefficient must be cancelled 
by the other coefficients in (3.14). We will see that to achieve this, two summations 
have to be done in analytic terms. 

The further factors that we find for fixed mi, m 2 , m 3 depend on 

Cl = 7T2,1+Vr 3j i = <7 2)3 - &2-7r2,2, C2 = Si-mi+7T 2j2 , CTjj = (Si-mi) + (Sj-mj) 

(5.17) 

and we have to evaluate (the replacement (a 2 ,pi) = —(02, P3) contributes to the 
phase factor) 

^ ' mi) \m 2 J \m 3 J \k 2 , 7r 2i2 , s 2 -m 2 -k 2 - tx %2 



^2,^2, 2 

(gi - mi + 7r 2i2 )!(q-2,3 - h- 7r 2;2 )! 
(A + 2- k 2 )\ 



>i, p 2 ) si - mi (a 2 , p 3 ) S2 ~ m2 (a 3 , Pi) S3 " m3 (5.18) 



where we must sum over k 2 and 7r 2j2 finally. The result is, after evaluating two finite 
2 Fi(l) series and applying the well known identity for T(z)T{l — z) 

1 , \ 1 / \ 1 / ^( s x \ ( s 2 \( s 3 



(A + 2)! 



si - mi)!( S2 - m 2 )\(s 3 - m 3 )! ( " l )( $2 )( * 3 ) Wen, p^)^ (5.19) 

\m 1 J \m 2 J \m 3 J AA 



Finally we take the remaining factors into account and cancel the three binomials 
in (5.11) after multiplication with the second factor of (5.16) with the result 



si!s 2 !s 3 ! ( s 3 
(2s 3 )!(A + 2)! V"i,2, "2,3, "1,3 



(5.20) 



which is the result expected from [7] up to the normalizing factor and coincides with 
(5.8). 
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5.1.2 The cases n > 1 

For N = 3 we use (3.4) and expand Z n where 2n is bounded by A 

ZH = E ( Al) T 2! A3 )-i Al+A3 -2 1+A2 r 3 A2+A3 (^) Al (^) A2 (^) A3 (5-21) 

Instead of (3.18) we obtain now after integration over the Tj and inserting M ra = 
M + 2n, A n = A — 2n, M = 5 — A = 2s 3 (we use primed quantities for m 8 and 
riij to mark the dependence on n) 

Mp + l ^ , ^ +n + n\ /M + 2n\ 1 



i,j partitions {tt g a },l 



l) n f n \. (Ai + A 3 )!(A 1 + A 2 + Ci)!(A 2 + A 3 + C2)! 



n\ VVAi,A 2 ,A 3 ; (S-M -K + 2)\ 

3 



1 2 

(-(a* K + «*_i))) fc < l[(a e ,w,r^ (5.22) 



7=1 



following the same arguments (5.3)-(5.5) as in the case n — 0. The group of terms 
we are considering is 

n<-i)*( m ^ 

= (-l)"^[n f Si ,V/ 2 V 01 ^ 2 )"' 2 ^ 2 '^)^ 1 ^ 2 ' 2 ^ 2 ^,^)^ 1 (5.23) 

Thus it remains to do the summations over k 2 and 7^,2- Using the results (5.3) - (5.7) 
we get first the polynomial in Wi (5.8) which is now of degree A n . The differential 
operator that is only a polynomial in the Laplacians results from the sum 

( -ir £ [(s - Mo -k 2 + 2)!]- 1 J2 ( ?r ( f r ( ?r 

: rx . rr, Ai! A 2 ! A 3 ! 

(Ax + A 3 )!(Ax + A 2 + Ci)!(A 2 + A 3 + C2)! (5.24) 

where (as in (5.17)) 

Ci = 0"2, 3 - k 2 - vr 2i 2, C2 = si - m[ + n 2:2 , = s i -m' i + Sj - (5.25) 
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The technique of summation is the same as in the case n = and we get altogether 



22 ( s i - m 'i + ^2 + A 3 )!(s 2 - m' 2 + Ai + A 3 )!(s 3 - m' 3 + Ai + A 2 )! 

{n( a -s) <9 }(«i^2) si - m ka 2 ,P3) s2 - m Ka3,Pi) s3 -^n^r (5 - 26) 

r,q 1=1 £ ' 



A certain simplification of this expression is possible by 

-1 3 

s 3 + zn \ 



x(s 2 - m' 2 + Ax + A 3 )!(s 3 - m' 3 + Ax + A 2 )! = {(2s 3 + 2n)\}- 1 s 1 \s 2 \s 3 \ 

x(si - m' + 1)a 2+ a 3 (s 2 -m' 2 + l) Al+A3 (s 3 - mf 3 + l)\ 1+ \ 2 (5.27) 

The sum over n of (5.26) vanishes, as long as divergences and traces are neglected, 
on shell and belongs to the kernel of the Noetherian, or can be eliminated by a field 
redefinition. The situation changes if we admit divergence terms. 

5.2 The divergence terms 

We allow for divergence terms 

(pd a )U s \p-a) (5.28) 

To extract trace field terms instead of divergence field terms would be more involved, 
since these arise from both the proper trace field and the deDonder field terms. In 
turn the deDonder field terms can be obtained by implementing a trace term into 
each divergence term as defined in (1.1), so that one trace field term results from 
each divergence term. But there are also the proper trace field terms. However, the 
divergence terms have also two different sources. 

5.2.1 The first source 

First we consider the group of terms from n = treated in subsection 5.1 and make 
the divergence terms explicit that were neglected there. In [7] the group of terms 
left over when we annihilate all traces, form the four blocks 



£(0,0) s £ (1,0) ; £(2,0) ^ £(3,0) {5 29) 
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traces=0 J 



Then equs. (2.2) and (3.16)-(3.18) of [7] can be reformulated as 

i=0 

i 

+ 5^(a i ,pi)(o i ,p i )- 



<9(a;,p m ) 
«9 2 



9(ai,pi +1 )9(a^,pj + i) 



+ (a 1 ,p 1 )(a 2 ,p 2 )(a 3 ,p 3 )- ^7 - w r}]£ ( °' 0) (5.30) 

<9(ai,p 2 )<9(a 2 ,p 3 )<9(a 3 ,pi) 

Comparison with (5.1) the curly bracket in (5.30) yields the differential operator 
part 

1 3 

2 5^ II[ n ^ ai '- Pi ^ <5i ( ai '- Pi + 1 ) n4 * (Pi+i l»=3=Pi)» n i = s i -m i (5.31) 

where <5j is either zero or one and Yli=% $i — 1- 

To obtain these divergence terms from the n = expression in subsection 5.1.1 
(see (2. 22), (5. 23)) we expand, keeping divergence terms 



K 2,Al 2,-^2 1 



^2/ V-^2,1/ \Al,2 



(0!, Pl )*l+^.l+Al, a ( a2j p2 ) A 2 ,l+«2 ( a3) p 3 )-3,2+fc3 

(a 2 ,p 3 ) 7r2 ' 1-A2 ' 1+7r2 ' 2+fe2 ~' e2 (a3,p0 7r3 ' 1 (ai,P2) 7ri ' 2 ~ Al ' 2 (5.32) 

We denote the powers of the divergence terms by 

5i = h + 71-1,1 + Ai j2 , 5 2 = k 2 + A 2) i, <5 3 = fc 3 + 7r 3i2 (5.33) 

Then the last line of (5.32) can, see (5.14), (5.15)), be written as 

(a2,p 3 ) S2 - m2 - S %a 3 ,p i y>- ms - s *(a 1 ,p 2 yi- m ^ (5.34) 

which agrees with (5.31). 

The operator form being correct, we set out to calculate the coefficients. We will 
do this here only for the case 

5i = 1,5 2 = 5 3 = 0. (5.35) 

If this coefficient comes out as predicted by (5.30), (5.31), we have no doubt that 
the coefficients in all other cases come out correctly, too. 
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We proceed as follows. From (5.33) we have only three choices for the parameters 

k ± e {0, 1}, TTx,! G {0, 1}, h ■ 7T lfl = (5.36) 

We treat these cases separately. Then in each case we have to do 2-^1 (1) summations 
in the way we did before. Let us formulate the program by 



{<5fcl,0^7Tl,l,0 
+ <Sfcl, 0^71-1,1,1 



(si - mi) 

fcl,fe2,fc3,7Tl,l,7r2,2,7r3,2 

o"2,3 -k 2 - vr 2 , 2 )!(si - m 1 + 7r 2)2 )! 



(A + 2-A; 2 )! 
o"2,3 - A; 2 - vr 2i2 + l)!(si - m 1 + 7r 2 , 2 - 1)! 



(A + 2-A; 2 )! 
1, , (q- 2 , 3 - k 2 - 7r 2|2 )!(si ~mi + tt 2 , 2 - 1)! 1 
2 (Ao + l-iba)! 1 

x(-l) m2+fc2 f, S2_m2 W°<W,o (5.37) 

with <7jj from (5.17). The summations are performed in the same fashion as in 
subsection 5.1 and the result is 



1 ( s 3 



) II( S< ~ 

^ i=i 



(2s 3 )!(A + 2)! Vm,2, ni >3 , n 2>3 

x^{(«i - m i) - ( s 2 - m 2 ) + (s 3 - m 3 )} (5.38) 
where the factorials cancel against factorials from 

11 (:,.) < 5 - 39 > 



=1 



Comparing this result with (5.20) we obtain as relevant factor the second line of 
(5.38). 

5.2.2 The second source 

This group of terms arises from n > 1. For an evaluation it is more practical to use 

(d ai ,d aj )(d avPj ) (5.40) 

acting on the fields h^(xi, Oj) as differential operators. We start from Section 5.1.1 
and our aim is to replace the Laplacians pf by the tracefree Fronsdal fields 

Ti = [~Pi + {a hPi ){ Ph d ai )]h^\p-a) (5.41) 
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In this fashion we must, after neglecting all Fronsdal fields (e.g. eliminating them 
by a field redefinition), replace the last factor in (5.26) by 



g(A)=n A = { A„A 2 ,A 3 } (5.42) 
1=1 e ' 

Because there is at most one divergence operation of each type (two of the same 
type vanish on shell and are, similarly as the Fronsdal operators, not resumed in the 
ansatz for a Lagrangian), of interest are only the cases 

Aj < 1 (5.43) 

selected from pole terms n up to three. 

Now we consider the case A = (1, 0, 0) which in (5.26) necessitates only to 
analyse the case n = 1. If we suppose that all Fronsdal operators are directly acting 
on the field h, it follows that the operators (di,Pi)(pi,d ai ) arising act also directly 
on these fields. This implies that the gradient operator, namely the first factor in 
this combination, stands to the right of all differential operators that appear in the 
same order as in (5.26) 

3 

r,q i=l 

Our strategy is then to move it to the left by commutators. At the left end we get 
zero because we get s + 1 derivations with respect to a on h^ s \ A commutator of 
the gradient (ai,pi) with (d ai ,p2) gives the sum of three Laplacians which lead to 
Fronsdal and further divergence fields. They can be neglected here. Thus we have 
only to deal with the commutator 

[Y[{da r ,d aq )<*,(a 1 ,p 1 )] = 

(d a2 , d a3 r's.3 k i2 (d ai , & az r'1.3 (d ai , d a2 ) n ^-\d a2 , Pl ) + 

n^(d ai ,d a J^-\d ai ,d a2 r'^(d a:i , Pl )} (5.45) 

Neglecting once again an additional divergence operator we replace (d a2 ,pi) by 
-(d a2 ,p3,). 

The next task is to identify the parameters {n^-,m^} by the {n^jjirii}. To do 
this we have to treat the two cases included in (5.45) separately. The term (I) (the 
second in (5.45)) is 

(/): M > s(9a 1 ,9« a ) Bi ' a (9a 1 ,9«.) Bi -»- 1 (9« a ,9«.)^" 

(d^ , p 2 ) S1 - m i (0 a2 , p,)^ {8 a , , Pi r~< +1 (5.46) 
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Comparison with (5.31), (5.34) yields 

n lj2 = n' 1)2 , ni >3 = n' 1)3 - 1, n 2 , 3 = n 2 3 ; 
mi = m^ — 1, m 2 = m 2 , m 3 = m 3 — 1; 
Si — mi — 1 — Si — m[, s 2 — m 2 = s 2 — m' 2 , s 3 — m 3 = s 3 — m 3 + 1 (5.47) 

The term (II) is 

(II) ■ ~n' h2 (d ai , d a2 f^- 1 (d ai , d a , )"i.3 («9 a2 , 9 a3 )"k 

(9 ai , P 2 ) si - mi (9 a2 , ps)*""** 1 (<9 a3 , pO*""* (5.48) 

Comparison with (5.31), (5.34) yields now 

"1,2 = rai.2 - 1, ni >3 = n' 1)3 , n 2j3 = n 23 ; 
mi = m[ — 1, m 2 = m 2 — 1, m 3 = m 3 
Si — mi — 1 = si — m' 1? s 2 — m 2 = s 2 — m' 2 + 1, s 3 — m 3 = s 3 — m 3 (5.49) 

All these relations are in both cases consistent in the sense 

Y n ^ = Y ~ lj Y mi = Y m 'i - 2 ( 5 - 5 °) 

i,j i,j i i 

The replacements 

Wij} ->• {Wij}, K} ->• {m,} (5.51) 
in the first line of (5.22) yield 

-M 



C}V (S3+1)C ™ 

(f n 2 ) " n f ''.) = i< 2s » + 2 » ! i" n < 5 - 52 > 

V / M v«v 7=i ( s j - m j)- 



3~- 

and for the whole first line in (5.26) 



lc^%[(2s 3 )\(A + 2)!]" 1 f[ (5-53) 



i=l 



From (5.27) we find that the summations lead to the factors 

(s 2 -m 2 + l)(s 3 -m 3 + l) (5.54) 

which for the case I equals 

(I): = (s 2 - m 2 + l)(s 3 - m 3 ) (5.55) 
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and for the case II is 

(//): = (s 2 -m 2 )(s 3 -m 3 + l) (5.56) 

These have to be subtracted from each other 

(/) - (//) : = -(s 2 - m 2 ) + (s 3 - m 3 ) (5.57) 

After multiplication with 1/2 we subtract this quantity from the last line in (5.38) 
and obtain 

-(si-mi) (5.58) 
thus getting the expected result \ri\ = |(si — mi) in (5.31). 

6 The quartic Noetherian 

In the case of N = 4 there are three inequivalent possibilities to contract the currents, 
which lead each to gauge invariant expressions. Only their sum with fixed special 
values for the three coupling constants exhibits the desired Bose symmetry. We 
denote the three different contraction schemes 

...Si <-> S 2 -H- S3 <H> S4 <H> Si... 
...Si <->■ S3 ■<->■ S 2 -H- S4 -H- Si... 

...Si f> S3 f> S 4 f> S 2 ^ Si... (6.1) 

by s-,t-, and -u-channel, respectively. We study only the s-channel expression, the 
other ones are obtained by corresponding permutations. For the leading terms we 
neglect once again the divergence terms 

(a t , Pt ) = (6.2) 

6.1 The case n = 

In this case n = we extract the factors 

so that from (3.14) remain as relevant factors 

y- (IIL1CJ) A, 1)V J s e -m e 

3 

(a^w^Hia^w^x* (6.4) 

a=l 
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By (6.2) this implies the following restrictions on the parameters fa, ii£ ta 

fa = fa = 0, 7T1,! = 7T 4 ,3 = (6.5) 

and the Fourier transformed differential operator part of (3.14) is 



1,3 



(a u w 2 y^ ( ai ,w 3 y^(a 2 , Wl y^ (a 2 , w 3 y^ 

(a 3 , Wl y^ (a 3 , w 2 y^ (a 4 , (a 4 , w 2 ) p ^ (6.6) 

where 

P2,l = 7T 2 ,1 + 7T 2i 2 + fc 2 , P3,2 = ^3,2 + ^3,3 + k 3 (6.7) 

whereas for all other cases 

Pij = ( 6 - 8 ) 

Moreover we note that 

Y,Pi,i =S-M = A, J2*i,J = S-M-K,K = k 2 + k 3 (6.9) 

Now we consider the p^j as fixed, leaving as variable parameters, say 

fa, fa, ^2,2, 7T3.3 (6.10) 

We have to sum over these. 

Next we express the quantities ( a by these now fixed variables p it j, respectively 
the variable parameters (6.10) 

Cl = P2,l + P3,l + P4,l - 7T2,2 ~ fa = R4 - 7T 2i 2 ~ fa 
C2 = Pl,2 + P3,2 + PA,2 + 7T2,2 ~ ^3,3 - ^3 = R-2 + ^2,2 ~ ^3,3 ~ fa 

C3 = Pl,3 + P2,3 + 7T 3 ,3 = R-3 + 7T 3 ,3 (6.11) 

We use the shorthands 

Si-mi = Vi (6.12) 

which allows us to express 

Y,U = A-k 2 -k 3 (6.13) 



and to present (6.4) in the form 

SlWs 3 W v / M \ 1 

partitions{mj } { PlJ } l ,0 ) / 1 li=l ■ 



V f 1 ^2+fc 3 ( R l~ 7r 2,2-fe2)!(R2 + 7r2,2-7r3,3-fe3)!(R3 + 7r3,3)! 

(A-fc 2 -fc 3 + 3)! 

^1 \ / ^2 \ / ^3 

Pl,2, Pl,3/ \fa, P2,l - 7r 2 , 2 - fa, 7T2,2, P2,3/ \^3, P3,l, P3,2 ~ ^3,3 ~ fa, 7C 3t3 



X 



\P4,1,P4,2/ " " 
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We recognize that the two quadrinomial coefficients depend on the free parameters 
n 2,2, 7T3,3, k 2 , k 3 . We denote this function within the curly bracket of (6.14) as 

^1,^2,^3,^4 la 

Pl,2,Pl,3,P2,l>P2,3,P3,l,P3,2,P4,l,P4,2 V / 

The summations over k 2 , k 3 , ^2,2, ^3,3 can be performed by reduction to Gaussian 
series 2-^1(1) for the first three summations and the fourth to a finite generalized 
hypergeo metric series 3-F 2 (l) of p3 )2 + 1 terms. The result is 

Pl,2,pl,3,P2,l,P2,3,P3,l,P3,2,P4,l,p4,2 
(_1)P2,1+P3,2 / Vx \ / V2 \ / U3 \ / ^ 



EL ^ ! \Pl,2, Pl,3/ \P2,1,P2,3J \P3,1, P3,2/ \P4,1, P4,2 

x (Rl " P2 (A + 3/ 3 " )!R3! (A " R2 + 3) <- (A " Ri " + 2 )^ 
x 3 i 71 2 (-P3,2, R3 + 1, A — Ri — R 2 + P2,i + 2; A - R x - R 2 + 2, A - R 2 + 3; 1) 

(6.16) 

Here we have not made use of the fact that Ri + R 2 + R3 = A, but we have kept 
these four parameters as independent. This turns out to be very practical when we 
treat the cases n > 0. 

6.2 The cases n > 1 

In this case we introduce the factor Z which in the notation using the momenta 
{wi} is 

3 

Z = J^r i+ i(l - r i+1 )w* -2^T i+1 T J+1 {w h w 3 ) (6.17) 

i=l i<j 

which can be expressed by the basis described in the Appendix 

4 

Z = ^A^Wi - Wi_i) 2 + B j;k (wj,w k ), (6.18) 

i=l l<j<fc<3 

where the coefficients Ai,Bj^ (8.1) are functions of T2,T3,t 4 specified in (8.5). 
Any n'th power of Z can be multinomially expanded 



ni,n2,{Xi},{nj,k} 



n \ / ni W n 2 
ni, n 2 / \Ai, A 2 , A 3 , A 4 / ^1,2, A*i,3, A*2,3 



(( Wl ) 2 ) Al ((^2 - ^l) 2 ) A2 (K - ^2) 2 ) A3 ((^3) 2 ) A4 

(w 1 ,w 2 r^(w 1 ,w 3 r^(w 2 ,w 3 r^A x ^A^A x 3 3 A^ B^ + ^B^ (6.19) 

The product over the A and 5 coefficients has to be integrated with the measure 

» 3 3 

(\[dT 3+1 e(T 3+1 ))\[rt +l (6.20) 
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where in analogy with (6. 3), (6. 4) we have 

4 3 

(a — 7Ti,a, Ca = A n ~ k 2 — k 3 , 

i=l o=l 

A n = A - 2n (6.21) 

and 

^ n^- = X -M n = l -M + n, ^ m' fc = M n = M + 2n (6.22) 

The terms n\ ^ belong to the kernel, whereas those with n\ — contribute to the 
leading terms. We will continue to discuss only these. 
In this case the number of derivatives A restricts n to 

< n < ^A (6.23) 

In such case we must integrate 

v ( n ) (w 1 ,w 2 r^(w 1 ,w 3 r^(w 2 ,w 3 r" i b^ + ^b^ = 

^ V/^1,2, ^1,3,^2,3/ 

J2( U ) (w 1 ,w 3 r^((w 1 ,w 2 ) + (w 2 , w 3 )r^ B?f + ^B^ (6.24) 
In the sequel we denote 

A*i,2 + A*2,3 = A*i,2,3 = n - A*i,3 (6-25) 

The integration is over the measure (6.20) and after introduction of the new 
variable T\ = 1 — t 2 — r 3 — r 4 yields (— 2) m - a times (see (3.5)) 



/ 



dTidT 2 aT 3 aT4 1 d{L — t± — r 2 — r 3 — r 4 jr 1 r 2 r 3 r 4 

= ^1,2,3K0 + /il,3)KC2 + /il,2,3)KC3 + ^gjj 

(A- fc 2 - A; 3 + 3)! 



(6.26) 



where the denominator is as usual independent of n. Taking into account the factor 
(— l) n /n\ from the residue of the gamma function (2.28) we obtain a finite sum 



(Cl+/il,3)!(C2 + ™-^ li3 )!(C 3 + /il,3)! 



,,,,;!( A - k 2 ~ /,, | 3)! 

(-2( Wl , W3 ))^ 3 [( Wl , W2 ) + (w 2 ,w 3 )] n -^ 3 (6.27) 

of the type of a (cut off at /ii j3 = n) Gaussian hypergeometric series. From now on 
we simplify the notations a little bit by defining 

A*i,3 = At, A*i,2,3 = n- fi (6.28) 
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For arbitrary n and using as before M n = M + 2n and corresponding labels 
{ni j,m' k } we obtain as starting point instead of (3.14) 

partitions {n^ .} Vl Vl kS / i<j 

v (c(+^)!(a+^-^)!(q+/i)! 

^ //!(A-fc 2 -fc 3 + 3)! 

partitions {n£ a },k( 

(l/2(a,, W£ + ^_ 1 ))^(J](a,,«; a )^)[-2( Wl , W 3)] M [K,^2) + (w 2 , w 3 )] n_/i (6.29) 

Here the parameters fc 2 , & 3 , 712,2, ^3,3 remain free and must be summed over. 

Now we proceed exactly as before and introduce a modified coefficient function 
C depending on n, p, v[ — — w! i 

f7"^2-">4 [„ ,,1 = \ (-l)k 2 +k 3 

P1,2,P1,3,P2,1,P2,3,P3,1>P3,2,/>4,1,P4,2L'"' h"! T-r4 / \ > 

1 li=l i k 2 ,k 3 ,^2,2,Tr3,3 

(Ri + p - 7r 2;2 - £: 2 )!(R 2 + n - p + tt 2i2 - 7r 3>3 - £: 3 )!(R 3 + P + 7r 3>3 )! 

//!(A - fc 2 - fc 3 + 3) 

\Pl,2, P\,3j \h, P2,l ~ ^2,2 ~ h, vr 2j2 , p 2 ,3/ \h, P3,l, P3,2 - 7r 3j3 - fc 3 , 7r 3i3 

xf *M (6.30) 

It is obvious that for general n, p we have only to replace in (6.14) and (6.16) 

Ri Ri + p, R 2 R 2 + n - p, R 3 =4> R 3 + p (6.31) 

whereas we keep A fixed, what is allowed if these four parameters are considered as 
independent. So instead of (6.16) we get now 



C u 'f "2 ^3 ."4 r n m1 — 

P1,2,P1,3,P2,1,P2,3,P3,1,P3,2P4,1,P4,2 L ' r^J 
(_1)P2,1+P3,2 / ^ \ / j/ 2 \ / ^ \ / ^ 



X 



^ ! \Pl,2, Pl,3/ \P2,1,P2,3/ \P3,1,P3,2/ \p4,l, 04,2, 

(Ri + /i - p 2 ,i)!(R 2 + n - // - P3, 2 )!(R 3 + p)! 



x 



p!(A + 3)! 

(A - R 2 - n + p + 3) P32 (A - Ri - R 2 - n + 2) P21 x 
3F 2 {-p3,2, R 3 + H + 1, A - Ri - R 2 - n + p 2 ,i + 2; 

A-R 1 -R 2 -n + 2,A-R 2 -n + /i + 3;l) (6.32) 

where < p < n < §. Thus for n = (6.32) reduces to (6.16). 
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Finally we give the result for all leading terms of the Lagrangian where the m' k 
are determined by the n.y through (2.44) 



A 




l,a i<j 



[-2(wi, U7 3 )] M (u>i + w 3 , w 2 ) 



(6.33) 



7 Conclusion 

For the N = 2 Noetherian we have obtained an expression which agreed with Frons- 
dal's free Lagrangian up to the trace terms, but the latter ones where wrong. The 
Z?-parameter had to be chosen zero in order not to produce .D-dependent coefficients 
of the trace terms. Since the N = 2 Noetherian is conformally invariant, the free 
Lagrangian cannot be. 

For the case iV = 3 we reproduced the leading terms and some of the divergence 
terms (deDonder terms for vanishing trace field) of the known cubic interaction 
Lagrangian [HEJE]. We know from [8] that these terms are sufficient to reproduce 
the whole interaction Lagrangian uniquely ("Proposition 1"), so that our present 
result agrees completely with our former result for the cubic interaction, and this 
interaction is now known to be conformally invariant. 

For the case N = 4 three versions (s-, t-, and u-channel) of the Noether potential 
exist. Bose symmetry requires a symmetric sum of them to be the complete expres- 
sion. We derived the leading terms of the s-channel version. It would be interesting 
whether these incomplete results can be completed by a version of "Proposition 1", 
perhaps provable from conformal invariance. 

The algorithm presented here is also applicable to analogous problems on spaces 
with a conformal structure such as AdS spaces. 

Acknowledgement: The author thanks Ruben Manvelyan and Karapet Mkrtchyan 
for clarifying discussions and helpful advices, and the organizers of the Third Work- 
shop on String Theory, July 2011, at the Witwatersrand University in Johannesburg, 
where a preliminary version of this work was presented. 

8 Appendix: A basis of bilinear forms for the mo- 
menta Wi 

In our analysis the Laplacians 



(wi - Wi-!) 2 , i G {1,2, ... 



N} 



(8.1) 
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play a special role for the selection of the kernel which vanishes on shell. Therefore 
we make the following ansatz for an arbitrary bilinear form 

N 

J2M^i-Wi-i) 2 + B i,i( w ii w i) ( 8 - 2 ) 

i=l l<i<j<N-l 

The dimension of this space of bilinear forms, namely N + (^J 1 ) is, however, one 
too big, since by momentum conservation we may set A N = 0. On the other hand 
we need all Laplacians. Therefore we try to imply a constraint on the coefficients 
Bij which is noncontradictory, e.g. the symmetric constraint 

-Bi,2 = B N -2,n-i (8.3) 

Applying this ansatz to Z (6.17), we obtain the equations 

A 1 + A 2 = t 2 (1-t 2 ) 
A 2 + A 3 = r 3 (l - r 3 ) 
A 3 + A 4 = r 4 (l - r 4 ) 
B lj3 = — 2r 2 r 4 
Bi >2 — 2A 2 = — 2r 2 r 3 

^2,3 - 2A 3 = -2r 3 r 4 (8.4) 
This leads to the following solutions 

Ai = r 2 (l - r 2 ) - ^r 3 (l + t 2 -t 3 - r 4 ) 
1 . 

A 2 = -r 3 {l + r 2 - r 3 - r 4 ) 
^3 = 7^3(1 - r 2 - r 3 + r 4 ) 

A 4 = r 4 (l - r 4 ) - ^r 3 (l - r 2 - r 3 + r 4 ) 

^i,2 = 5 2i3 = r 3 (l - r 2 - t 3 - r 4 ) 

fii,3 = -2r 2 r 4 (8.5) 
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